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Drifting pattern domains (DPDs), i.e. moving localized patches of traveling waves embedded in 
a stationary (Turing) pattern background and vice versa, are observed in simulations of a reaction- 
diffusion model with nonlocal coupling. Within this model, a region of bistability between Turing 
patterns and traveling waves arises from a codimension-2 Turing- wave bifurcation (TWB). DPDs are 
found within that region in a substantial distance from the TWB. We investigated the dynamics of 
single interfaces between Turing and wave patterns. It is found that DPDs exist due to a locking of 
the interface velocities, which is imposed by the absence of space-time defects near these interfaces. 

PACS numbers: 03.40.Kf, 47.54. +r, 82.40.Ck 
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- Introduction - Pattern forming processes in nonequi- 
librium systems can be classified according to the pri- 
mary instability of the spatially homogeneous state. Ref. 
jjj distinguishes three basic types of instabilities in un- 
bounded systems: (i) spatially periodic and stationary 
in time, (ii) spatially periodic and oscillatory in time 
and (in) spatially homogeneous and oscillatory in time. 
Within the reaction-diffusion literature, these instabil- 
ities are known as Turing, wave and Hopf bifurcation, 
respectively. 

Many chemical and biological patterns are well cap- 
tured by so called activator-inhibitor models || describ- 
ing the dynamics of two reacting and diffusing substances 
with two coupled partial differential equations. In such 
two component reaction-diffusion models only Turing 
and Hopf instabilities are possible. Recently, numeri- 
cal investigations of chemical reaction-diffusion systems 
with three components j^] and nonlocal coupling Q have 
yielded the occurrence of wave instabilities and the cor- 
responding patterns. A universal description of patterns 
near these instabilities is achieved within the framework 
of amplitude equations . 

Here, we study a simple FitzHugh-Nagumo model with 
inhibitory nonlocal coupling that is obtained as a limiting 
case of a three component reaction-diffusion system. It 
describes the interaction of an activator species with an 
inhibitor. For slow inhibitor diffusion (compared to the 
activator diffusion), the model exhibits wave instabilities, 
while, for fast inhibitor diffusion, Turing instabilities are 
found. The two instabilities occur simultaneously at a 
codimcnsion-2 Turing- wave bifurcation (TWB). Such a 
situation has been found earlier within a model for bi- 
nary convection || and is a generalization of the well 
investigated Turing-Hopf instability in reaction-diffusion 
systems jfj . Basic properties of a TWB have been studied 
theoretically in amplitude equations Q as well as exper- 
imentally in a one-dimensional gas-discharge system 
In our model, we find a pattern previously unknown in 




FIG. 1. Space-time plots of the field u in grey scale for 
three examples of DPDs found in numerical simulations of 
Eqs. (1). In all three cases a — 6.0, and only a part of the 
system of length L = 409.6 is shown, (a) Large DPD with 
S = 0.84. (b),(c) DPDs consisting of a single cell of Turing 
and wave respectively. In (b) S — 0.80 and in (c) 5 = 0.91. 
Other parameters, see 11 7f| - 



reaction-diffusion systems: drifting pattern domains 
(DPDs), i.e. localized patches of traveling waves em- 
bedded in a Turing background and vice versa (see Fig. 
1). These patches have constant width and move (drift) 
with constant speed. As they drift along, maxima of the 
concentrations of activator and inhibitor are conserved 
beyond both boundaries of the DPD, where Turing and 
wave patterns are joined together. Consequently, for- 
mation of space-time defects by coalescence of maxima 
and minima is prevented. Similar patterns have been 
reported in a variety of hydrodynamical experimental 
systems, see e.g. |l(]|n]] and have been related to sec- 
ondary instabilities (parity breaking) of stationary pat- 
terns DPDs exist in a broad region of the para- 
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FIG. 2. (a) Critical wavenumber against inverse non- 
local coupling range a. (b) Real (x) and imaginary part (w) 
of X(k) at the TWB, for parameters see Q. 

meter space, but appear only in a substantial distance 
to the onset of pattern formation. Their existence re- 
gion is characterized by bistability between waves and 
Turing patterns. Near the boundary of DPD existence 
small DPDs containing only a single wave or Turing 
„cell" inside a background of the respective other state 
are encountered (see Figs. lb,c). Outside the DPD ex- 
istence region, the width of pattern domains shrinks or 
expands; these transient domains exhibit defects at the 
interfaces. Large DPDs are composed of two interfaces 
separating wave and Turing patterns. In this Letter we 
will study the dynamics of such interfaces. These in- 
terfaces typically select the wavenumber of the invading 
domain and form one parameter families characterized 
by the wavenumber of the invaded domain. Far away 
from the TWB, interfaces can exhibit a locking mech- 
anism of their velocities due to the absence of defects. 
This locking implies that the interface velocity is fixed 
by the wavenumbers and frequency of the patterns that 
they separate. If both interfaces in a DPD are locked, 
they have to travel with equal speed. This mechanism 
is responsible for the existence of DPDs with constant 
width. 

- Model equations and linear stability - We study a 
variant of the FitzHugh-Nagumo equation supplemented 
by an inhibitory nonlocal coupling in the dynamics of the 
activator u 



dtu = au + f3u 2 — au 3 — bv + d 2 u 

-a\x—x | 

dtv = cu — dv + Sd x v. 



H I e " ~ 'u(x , t)dx 



(1) 



Eqs. (1) represent a limiting case of a three variable 
model involving the activator u, the inhibitor v and an 
additional fast inhibitor |l3| . Related three variable mod- 
els have been introduced previously to describe pattern 
formation on sea shells and in cell biology H] as well 
as spot dynamics in gas discharges Jl5[ and concentra- 
tion patterns in heterogeneous catalysis |h|. Here, the 
emphasis is on the onset of pattern formation result- 
ing from destabilization of a single homogeneous steady 
state. Eqs. (1) possess the trivial homogeneous fixed 



point u = (uo, vq) t — (0, 0) T for all parameter values. 
Here, we consider the regime where this fixed point is 
the only one present, i.e. a < bc/d + 2/x/cr and consider 
perturbations oc e lkx ~ •M fc )* 5 where X(k) = x(fc) + iuj(k). 
The growth rates X(k) are given by the eigenvalues of 
the Jacobian. Linear stability analysis reveals that Eqs. 
(1) exhibit wave instabilities if the nonlocal coupling is 
of sufficiently long range a < a c = (2/i/(l + S)) 1 ^ 3 . 

In the following we vary the control parameters a and 
6; the ,, driving force" a represents the kinetics, whereas 
the ratio of diffusion coefficients S describes the spa- 
tial coupling in the medium. All other parameters of 
Eqs. (1) have been fixed llTj]. For the wave bifurcation, 
the critical wavenumber k w and parameters aw , $w are 
obtained from the condition X(k w ) = ±iujQ where the 
perturbation with k w is the fastest growing mode with 

{k w ) 2 = yff| - o- 2 . Note, that for both a = a c and 

(7 = (global coupling limit) the critical wavenumber 
is k w = 0, see Fig. 2a. Similarly, a competing Tur- 
ing instability appears for a critical parameter with a 
wavenumber k^, where the leading eigenvalue A(fey) = 0. 
For large enough driving a, the wave instability appears 
for small <5, while for large S the Turing instability desta- 
bilizes the homogeneous state. For the chosen parameter 
values, the system exhibits a TWB point (see Fig. 3a 
and p7|]). For the corresponding A(fe), see Fig. 2b. 

- Weakly nonlinear analysis - Near the TWB, we can 
write u=(u, v) T as a perturbative expansion around u 
using a small parameter e, indicating the distance to the 
instability threshold: u = u + eui + e 2 u 2 + e 3 u 3 + • • • 
and use the following multiple scale ansatz: Q 

ui = [A(X,T v T 2 )\J A e i ( LUot+k '«' x) + B(X,T v T 2 )U B e l(ul0t - k ™ x) 
+K{X,T v T 2 )Yke lk T* + c.c.}/2. 

This leads to a set of coupled equations for the ampli- 
tudes A 7 B and 1Z for left-, right- going waves and Turing 
pattern that depend on slow time and space variables. 
After reestablishing the original time and space variables 
and performing further e-independent scaling, one ob- 
tains: 

d t Tl = t]11- \TZ\ 2 TZ + id 2 x Tl - c(|^i| 2 + \B^)Tl 
d t A + c g d x A = P A + (1 + ic x )d 2 x A - (1 - ic 3 )\A\ 2 A 

-g(l - ic 2 )\B\ 2 A - v{\ - in)\Tl\ 2 A 
d t B - c g d x B = P B + (l + i Cl )d 2 x B - (1 - ic 3 )\B\ 2 B 

-g{l-ic 2 )\A\ 2 B -is{l-iK)\1l\ 2 B. (2) 

For the detailed values of all coefficients, see Note, 
that the nonlocal term of Eqs. (1) only enters into the 
diffusion coefficients of Eqs. (2) and does not give rise 
to a nonlocal term in Eqs. (2). Knowledge of the coef- 
ficients of Eqs. (2) allows analytical predictions of the 
pattern dynamics. Here, traveling waves are always pre- 
ferred over standing waves (g > 1, see pj) and bistability 
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between wave and Turing patterns is found (vQ > 1). In 
this bistability region in parameter space (see Fig. 3a), a 
family of stable Turing patterns and two families of stable 
left- and right-traveling waves parametrised by their cor- 
responding wavenumbers coexist. To get further insight, 
we take a closer look at the dynamics of single interfaces 
separating domains of Turing and wave patterns. 

- Interface dynamics. - With suitable initial condi- 
tions, a moving interface between Turing and wave pat- 
terns will be formed in simulations of Eqs. (1). We can 
distinguish two types of interfaces depending on whether 
the phase velocity of the waves points towards the inter- 
face or away from it. This classification is independent 
from the direction in which the interface is moving. In 
the following, we will call the first type inward-interfaces 
and the second outward-interfaces. Fig. 3b and 3c show 
examples of the latter type. 

Near the TWB, we have studied general properties 
of such interfaces in amplitude equations (3) by count- 
ing arguments 0] as well as by direct numerical simula- 
tions. Counting arguments are applied to ordinary dif- 
ferential equations obtained from a coherent structure 
ansatz in a comoving frame. We observe that, typi- 
cally, the wavenumber of the invaded domain remains 
constant, while it adapts in the invading domain. In 
other words, the interface selects a particular wavenum- 
ber for the invading state, while the initial wavenumber 
of the invaded state is a free parameter. The velocity of 
the interface is a function of this parameter. For Turing 
patterns the selected wavenumber is always the critical 
one, i. e. k^ 1 — kj*, while for waves typically kyy ^ k^, 
and therefore ui sel ^ luq. This is valid for both inward- 
and outward-interfaces. Thus, we typically have two one- 
parameter families of interfaces for a given point in pa- 
rameter space. These results are confirmed by numerical 
simulations of the nonlocal model (1) near the TWB. 

Simulations of interfaces in Eqs. (1) far away from the 
TWB, show qualitatively similar behavior with respect 
to the selected wavenumbers. In addition, interfaces far 
away from the TWB may exhibit a locking mechanism of 
their velocities, which are fixed by the wavenumbers and 
frequencies of the Turing and wave domains. More specif- 
ically, the selected velocity is determined by the absence 
of defects at the interface. For geometrical reasons an in- 
terface without defects, that connects a wave state with 
wavenumber kyy and frequency u> and a Turing state with 
kr, has a speed |n; c/c| = (^/(kr—kw)- This velocity lock- 
ing mechanism is found for both types of interfaces. Two 
examples of areas where locking occurs (locking tongues) 
in parameter space for inward- and outward-interfaces 
are given in Fig. 3a. A locked outward-interface is dis- 
played in Fig. 3b. Outside the tongue the outward- 
interfaces display phase slips (see Fig. 3c). The area 
of the locking tongues depends only weakly on the front 
parameter kw for inward-interfaces and kx for outward- 
interfaces. Note, that the locking tongues open at a sub- 
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FIG. 3. (a) Parameter space 
(black circle). The light grey region indicates bistability be- 
tween the Turing and wave pattern with k^, k^y as predicted 
from Eqs. (2). Dark grey regions correspond to two exam- 
ples of locking tongues for an outward-interface with kr = k^ 
(region OT) and an inward-interface with kw = k^ (region 
IT). These two tongues are shown only up to a ~ 5.6. The 
dashed lines show where the selected velocities of interfaces 
in Eqs. (2) coincide with i>( oc fc. White circles indicate pa- 
rameter values of simulations shown in (b) and (c). See Fig. 
4 for a description of the dashed regions, (b) and (c) show 
space-time plots of u in grey scale from simulations of Eqs. 
(1) showing outward-interfaces for a = 5.2. In (b) an exam- 
ple inside the locking tongue is shown for 8 = 0.80 and in (c) 
an interface outside the tongue exhibiting defects (inside the 
white circles) for 8 = 0.86 is shown. 



stantial distance from the TWB. Since the rapidly vary- 
ing space and time scales have been factored out in the 
amplitude equations (2), locking tongues cannot be found 
therein. However on a line in the parameter space (see 
Fig. 3a), the velocity of interfaces in Eqs. (2) coincides 
with the velocity prescribed by the locking mechanism. 
The locking mechanism arises when the characteristic 
width of the interfaces is of the same order than the char- 
acteristic length scale of the patterns. 

- DPDs and their Phase Diagram. - Consider that 
a large DPD is composed of an inward-interface and an 
outward-interface, that practically do not interact (see 
e.g. Fig. la). If both interfaces exhibit no defects and 
are locked, their velocities have equal magnitude |uiocfc| 
but opposite signs. This ensures constant width and al- 
lows construction of DPDs of arbitrary size. Indeed, the 
region of existence of large DPDs starts to open where the 
locking tongues for both interface types begin to overlap 
(see Fig. 3a). This is the case for a > 5.7. Above that 
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FIG. 4. The region of existence of DPDs obtained in sim- 
ulations of Eq. (1) in the a-5 parameter space shown dashed. 
The light grey area corresponds to the bistable region between 
the critical wave and Turing patterns as calculated from the 
amplitude equations (3). In region B DPDs of any size exist; 
the size being determined only by the initial condition. In 
region C small domains of wave patches traveling in a Tur- 
ing background are found. In region A only Turing-droplets 
are stable. The three circles correspond to the location of 
simulations shown in Fig. 1. 

value, DPDs spontaneously form from a variety of ini- 
tial conditions. We have determined the parameter re- 
gion, where they propagate with constant width and drift 
speed, from extensive simulations in systems with sizes 
L > 400 and periodic boundary conditions. The results 
are shown in the phase diagram of Fig. 4. 

We can distinguish three different subregions. In re- 
gion B, DPDs of any size, with two locked interfaces 
traveling at the same speed, are found (see Fig. la). 
In region A, the inward-interface is no longer locked and 
its speed is smaller than |i>i oc /c|- Therefore large domains 
of Turing (wave) patterns contract (expand) in size until 
only a stable DPDs containing a single Turing cell is left 
(see Fig. lb). In region C, the outward-interface selects a 
kyy which would be unstable against Turing patterns in 
an infinite domain. Therefore, the wave domain forming 
the DPD is mostly replaced by a Turing pattern. How- 
ever, small DPDs with a few wavelength of wave pattern 
are still encountered. At the outer boundary of region C, 
only DPDs with a single wave cell are found to be stable 
(see Fig. lc). 

- Conclusion - We found a large variety of drifting pat- 
tern domains in a reaction-diffusion model with nonlocal 
coupling. Their ingredients include a bistability between 
wave and Turing patterns near a codimension-2 point as 
well as absence of defects at the interface. They exist 
as robust patterns only in a finite distance to the onset 
of pattern formation. Our results are not limited to the 
reaction-diffusion model studied here and should carry 
over to other physical systems with similar pattern form- 
ing instabilities. Altogether, DPDs and their constituting 
interfaces represent a generalization of simpler structures 
such as fronts and pulses in bistable reaction-diffusion 



systems, which do not simply combine two homogeneous 
states, but, instead, select their constituents from whole 
families of possible traveling or stationary periodic pat- 
terns. 
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